Chirality-driven mass enhancement in the kagome Hubbard model 
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We investigate the quasiparticle-mass enhancement in the Hubbard model on the frustrated 
kagome lattice by using a cluster extension of the dynamical mean-held theory. By analyzing the 
cluster density matrix, we Hnd a hierarchy of energy scale among charge, spin, and chirality degrees 
of freedom. Large amount of entropy associated with the chirality is released at a much lower tem- 
perature than other energy scales for spin and charge fluctuations, leading to a sharp peak in the 
speciflc heat and the single-particle spectrum. The results manifest a generic mechanism of mass 
enhancement driven by an emergent composite degree of freedom under geometrical frustration. 
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Heavy-fermion behavior is one of the most intriguing 
phenomena in correlated electron systems. Canonical ex- 
amples are found in rare-earth compounds, in which a 
large quasiparticle-mass enhancement is observed, e.g., 
in the specific-heat coefficient 7. The origin has been ar- 
gued on the basis of the Kondo effect — the screening of 
local /-electron moments by conduction electrons [Ij] . In 
this case, the local /-electron moments serve as an "en- 
tropy reservoir" for the heavy- mass behavior 0. 

Recently, a class of transition metal compounds has 
drawn considerable attention due to their heavy-fermion 
behavior. A large 7 of the order of lOOmJ/molK^ 
was observed in several compounds such as LiV204 :_3], 
Y(Sc)Mn2 [il, and /3-Mn The conventional Kondo 
scenario does not apply straightforwardly to this heavy- 
fermion behavior, since these Sd-electron compounds do 
not have an obvious entropy reservoir like the localized 
/ moments in the rare-earth compounds. There have 
been maiiy proposals for the origin of this puzzling be- 
havior |6l-[ll|. but the issue remains controversial so far. 

One of the proposals is the effect of strong electron cor- 
relation under geometrical frustration of the lattice struc- 
ture 0, S ■ Divergence of the quasiparticle mass at 
the correlation-driven Mott transition was firstly pointed 
out by Brinkman and Rice (l2j . and later, a more sophis- 
ticated picture was provided by the dynamical mean- 
field theory (DMFT) [3- It was shown that the elec- 
tron correlation results in separation of energy scale be- 
tween charge and spin (and orbital) degrees of freedom 
of electrons: Namely, it suppresses charge fiuctuations 
at an energy scale of the Coulomb repulsion, and leaves 
large spin (and orbital) fluctuations at much lower tem- 
peratures T, which serve as an entropy reservoir jT3 |. 
This mass divergence, however, persists only when spa- 
tial correlations are neglected. Indeed, when DMFT 
is extended to include spatial correlations, an antifer- 
romagnetic spin correlation develops at low T, which 
suppresses the local spin fiuctuations and collapses the 
heavy-fermion state 15| . It is widely considered that this 
obstacle can be circumvented by the geometrical frustra- 
tion: The frustration suppresses the development of the 
spatial correlation and rejuvenates the masked heavy- 
fermion behavior. Hence, in this scenario, the heavy- 



fermion behavior is an intrinsic property associated with 
the criticality of the Mott transition, and the geometri- 
cal frustration plays a secondary role of uncovering it by 
suppressing spatial correlations. 

In this Letter, contrary to the prevailing view, we re- 
veal an intensive role of the geometrical frustration on 
the quasiparticle-mass enhancement. We show that the 
frustration not only suppresses the spatial correlations 
but also brings into being a composite degree of free- 
dom, which plays a role of the entropy reservoir. It is 
known that a composite object, such as spin chirality jl^ 
and self-organized clusters often emerges in in- 

sulating systems, but we here explore an importance of 
such composites in correlated metals. By studying the 
kagome-lattice Hubbard model by a cluster extension of 
DMFT, we demonstrate that the spin chirality degree 
of freedom emerges even in the metallic region. We find 
that there is a hierarchy of energy scale among the charge, 
spin, and chirality, and that the chirality is dominant at 
the lowest T with carrying a large amount of entropy. 
The release of this chirality-associated entropy leads to 
the mass enhancement, with yielding a sharp peak in the 
specific heat as well as in the single-particle spectrum. 

We consider the Hubbard model on the kagome lattice 
shown in Fig. [TJ whose Hamiltonian is given by 



(1) 



where the sum of (ij) in the hopping term is taken over 
the nearest-neighbor sites, and riia = cj^Cia in the onsite 
Coulomb repulsion term. The model has been extensively 
studied as a minimal model including both electron cor- 
relation and geometrical frustration [2l| - |23 |. Hereafter 
we set t = I and the Boltzmann constant fee = 1, and fix 
the electron density at half filling. 

We adopt the cellular dynamical mean-field theory 
(CDMFT) [lOl to study this model, which is an ex- 
tension of DMFT to include the effect of spatial corre- 
lations within a cluster. The method is basically the 
same as that used in Ref. [23| . however, in the present 
study, we perform the CDMFT calculations for two dif- 
ferent clusters with three and nine sites (Fig. [T]) to check 
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the cluster-size dependence. Moreover, as a solver for 
the effective cluster model, we employ the continuous- 
time auxiliary-field quantum Monte Carlo (MC) method 
based on the perturbation expansion in terms of U [25| . 
This solver gives precise results more efficiently com- 
pared with the conventional Hirsch-Fye algorithm used 
in Ref. [1^. Typically we take 256 imaginary-time slices 
and 10^ MC steps, and the convergence is reached after 
20 self-consistency loops, e.g., at C/ = 6 and T = 0.1. 
Larger computational cost is required at lower T and 
larger U, in particular, in the calculations of the specific 
heat and the entropy (see below). The minus sign prob- 
lem is not severe in the parameter range that we have 
investigated. 

In order to identify the relevant degrees of freedom in 
the system, we calculate the cluster density matrix, which 
gives the probability distribution of quantum mechanical 
states within the cluster. Here we consider the density 
matrix defined on the three sites within each cluster in- 
dicated in Fig. [T] We calculate the diagonal components 

p* = lTr(|*)(vI/|e-«/^) (2) 

for a cluster state | ^t) , where Z is the partition function. 
The 43 = 64 states are classified by irreducible represen- 
tations under the U{1) (8) SU{2) (g) C^y symmetry of the 
model. Among them, we focus on pq, for the following 
three states (and their symmetrically equivalent ones); 

\D) = -i=(in,t,o) + |o,n,t) + it,o,n», (3) 

\S) = |t,t,t>, (4) 

\K) = -i=(|;,t,t) + '^lt4,t)+'^'lt,t,;)), (5) 

where | ti, ti 0) = cJ-^c|j^C2||vac), etc. (|vac) is a vacuum), 
and w is a phase factor related with the helicity (defined 
below). Here, we call \D) the doublon state which in- 
cludes a doublon- holon pair, \S) the spin-polarized state 
with the total spin S — 3/2, and |A') the chiral state. 
Among these states, the chiral state is of special interest. 
I A') retains fourfold degeneracy with its time- reversal and 




FIG. 1: (Color online) (a) Schematic picture of the mapping 
of the kagome lattice to three and nine-sites clusters in the 
CDMFT calculations. The three-site density matrix is cal- 
culated for the numbered three sites in each cluster, (b) 
Schematic pictures of the doublon, spin-polarized, and chi- 
ral states defined in Eqs. ([3])-([5]). 



reflection conjugate states: The four states are labeled by 
the two discrete quantum numbers, the z component of 
total spin 5'°* = ±1/2 and the helicity 77 = ±1 [tj in 
Eq. ([S]) is given by exp(i|7r77)]. We note that the chiral 
states are the eigenstates of scalar- and vector-chirality 
operators, studied in the previous CDMFT study p8l |. 

Figure [2] shows our CDMFT results for the density ma- 
trices, pd, PSi and px, in correlated metallic region at 
U = 6. (The critical value for the Mott transition was 
estimated at Uc — 8.2 ^^.) The three components show 
qualitatively different T dependences. At high T enough, 
PD ~ ps ~ PK ~ 1/64, since all the 64 states have nearly 
equal weights. As lowering T, pjy is suppressed but ps 
and pk are enhanced, ps turns to decrease with showing 
a broad peak at T ~ 1.2. px continues to increase down 
to lower T, but finally turns to decrease with showing a 
peak at r ~ 0.32. It is noteworthy that the results show 
little cluster-size dependence as shown in Fig. [5J This in- 
dicates that our CDMFT results show good convergence 
to the thermodynamic limit. Such rapid convergence may 
be attributed to suppressed inter-triangular correlations 
under the corner-sharing topology of the kagome lattice, 
as discussed in the localized spin models [3 [13] . 

These characteristic T dependences indicate a hier- 
archy of energy scale for relevant degrees of freedom. 
Charge fluctuations are firstly suppressed by large U, be- 
ing signaled by the suppression of which includes a 
doublon-holon pair. At a lower energy scale where ps 
decreases steeply, the spin degree of freedom is frozen 
out, and finally, the chirality is quenched at the low- 
est energy scale, corresponding to the suppression of pK. 
We identify the characteristic energy scales for charge, 
spin, and chirality degrees of freedom by the tempera- 
tures Tchargo, Tspin, and Tchirai where po, Ps, and pk are 
suppressed most rapidly (the inflection points indicated 
by the triangles in Fig. ^ . The hierarchy is clearly seen 

as Tchargo — 2.0 > Tspin ~ 0.37 > Tchirai — 0.18. 

We investigate the energy hierarchy systematically in 
the correlated metallic region. Figures [31^a)-(c) show T 




FIG. 2: (Color online) Temperature dependences of the three- 
site density matrices, po, ps, and pk- The results for the 
three(nine)-site cluster are shown by circles (crosses). 
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FIG. 3: (Color online) Temperature dependences of the den- 
sity matrices (a) pD, (b) ps, and (c) pK- The results are for 
the three-site cluster and the lines are guides for the eye. The 
characteristic energy scales Tchargc, rapin, and Tchirai, deter- 
mined by the steep decrease of po, ps, and pK, are plotted 
as a function of U in (d). The dashed line shows the effective 
exchange interaction J — At^ /U for comparison. 



dependences of po, Ps, and pK with varying U . The es- 
timated characteristic temperatures are summarized in 
the crossover phase diagram in Fig. [Hd). Tchargc in- 
creases as [/, which is reasonable since charge fluctu- 
ations are suppressed below T ~ t/ — W* {W* is a 
renormalized band width). On the other hand, Tspin de- 
creases as U increases. We compare Tgpin with an ef- 
fective superexchange interaction J = At'^/U, derived in 
the strong-coupling expansion in t/U, and find approx- 
imately Tgpin ~ J/2. This suggests that the reduction 
of Ps is a strong-correlation effect, originating from the 
effective spin exchange under suppressed charge fiuctua- 
tions. Tchirai also decreases as U increases, and is always 
located well below Tgpin. This indicates that the com- 
posite chiral degree of freedom is dominated by a much 
smaller energy scale compared with J. 

The energy hierarchy and the quenching of each degree 
of freedom are also observed from the entropic point of 
view. Figures Hlja) and IH^b) show the CDMFT results 
for the specific heat C and the entropy S, respectively. 
Here, we calculated C by numerical differentiation of the 
internal energy: To obtain the high precision data, we 
took a large number of imaginary-time slices, e.g., 8192 
shces at T ~ 0.1 and L/ ~ 6 dil. The entropy S is ob- 
tained by the numerical integration of C/T. As a result, 
C shows a sharp peak at low T ~ Tchirai in addition to 
a broad peak at T ^ Tchargc- The broad peak comes 
from the release of charge entropy; the peak temperature 
shifts to higher in accord with Tcharge as increasing U 
(not shown). On the other hand, the low-T sharp peak 
originates from the entropy release associated with the 




FIG. 4: (Color online) Temperature dependences of (a) spe- 
cific heat and (b) entropy a.t U — 6 per site. The results at 
[7 = are shown for comparison. Three vertical dashed lines 
indicate the characteristic temperatures estimated in Fig. [S] 
and three horizontal lines in (b) show the values of log 4/3, 
log 8/3, and log 64/3. (c) Single particle spectrum at i7 = 6 
obtained by the maximum entropy method, (d) Tempera- 
ture dependence of the uniform magnetic susceptibility, (e) 
Enhancement factors of the specific-heat coefficient estimated 
from Tchirai, 7*, and of the uniform magnetic susceptibility at 
T = 0.05, X* ■ See the text for details. All the results are 
obtained by the three-site cluster DMFT. 



chirality: In fact, as shown in Fig. lU^b), the entropy S 
takes a value close to log 4/3 at Tchirai, corresponding to 
the fourfold degeneracy of the chiral states \K). It is dif- 
ficult to see a clear structure in C at T ~ Tgpin, however, 
S reaches nearly log 8/3 at Tgpin, reflecting the fourfold 
spin-polarized states \S) in addition to the four chiral 
states. Our results indicate that the entropy associated 
with chirality, spin, and charge degrees of freedom is re- 
leased at separated energy scales Tchirai, Tspin, and Tchargc 
in Fig. [3Jd). To our knowledge, the resultant sharp peak 
in the specific heat due to the chirality has never been 
found in the previous studies. 

The sharp peak is considered as a hallmark of the for- 
mation of Fermi liquid state below Tchirai- In fact, as 
shown in Fig.[?I^c), the single-particle spectrum develops 
a sharp coherence peak at the Fermi level a; = below 
Tchirai- In addition, as plotted in Fig. IHJd), the uniform 
magnetic susceptibility x approaches a nonzero value as 
T — > after showing a broad peak at around Tchirai. 
Here x is calculated from the magnetization by applying 
a small magnetic field. All the results consistently indi- 
cate that a Fermi liquid state emerges at low T and the 
effective Fermi temperature is set by Tchirai - 

In the low-T Fermi liquid state, the quasiparticle mass 
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is largely enhanced because the effective Fermi tempera- 
ture Tchirai is suppressed. Although it is difficult to esti- 
mate the specific-heat coefficient 7 directly from the low- 
T data of C, we obtain a rough estimate of 7 from the 
value of Tchirai: Assuming that the chiral entropy log 4/3 
is released uniformly in the range of < T < Tchirai, 
we obtain the mass enhancement factor 7* = 7/70 — 
log4/(3Tchirai7o) — 3 at C/ = 6, where 70 — 0.873 is the 
specific-heat coefficient at [/ = 0. As Tchirai is suppressed 
for larger U [Fig. [2Id)] , 7* is enhanced with increasing 
U , as shown in Fig. SJe). We also plot the enhancement 
factor X* = x/xo estimated at T = 0.05, where xo is the 
susceptibility at [7 = 0. The results of 7* and x* consis- 
tently indicate a chirality-driven mass enhancement with 
approaching the Mott transition. 

In the formation of the heavy-quasiparticle state, the 
geometrical frustration plays an essential role, not only 
to suppress antiferromagnetic correlations but to bring 
about the degeneracy due to high local symmetry — in 
the present model, the fourfold degeneracy of the chiral 
states in the triangular unit. In general, geometrically- 
frustrated lattices possess such local unit with high sym- 
metry, e.g., triangles or tetrahedra, and hence, it is 
commonly expected that similar locally-degenerate state 
emerges. Our finding provides a prototypical example 
for the ubiquitous phenomenon in frustrated systems un- 
der strong correlation, that is, crossover from highly- 
symmetric local states to Fermi liquid state with an en- 
hancement of quasiparticle mass. We believe that the 
local degeneracy and emergent composite degree of free- 
dom are of primary importance for the puzzling heavy- 
fermion state in transition metal compounds, such as 
LiV204 3, Y(Sc)Mn2 [|, and /3-Mn Q. 

Finally let us remark on the fate of Tchirai as the system 
enters Mott insulating region at larger U . In the sense 
that it sets the effective Fermi temperature, Tchirai might 
terminate at the first-order Mott phase boundary. We 
note in fact that Tchirai in Fig. [31[d) appear to approach 
the critical point {Uc, Tc) suggested in Ref. [23|. On 
the other hand, the spin chirality itself can be important 
even in the insulating region: Higher order perturbations 
in terms oi t/U give higher-order exchange interactions 
in the effective spin model, which may induce compos- 



ite spin objects, such as the chirality. Unfortunately it 
is hard to clarify how Tchirai behaves at C/ > 8 because 
of the limitation of our simulation, and some comple- 
mentary study is highly desired. An interesting issue 
related to this problem is the recent experimental ob- 
servation of a sharp peak in the low-T specific heat of 
■^He on graphite at a commensurate filling 4/7 [33|. It 
has been argued that the system is insulating but in the 
vicinity of the Mott transition [U, and the geometry 
of "^He atoms has frustration in between triangular and 
kagome lattices [32:]. 

In summary, wc have studied the correlated metallic re- 
gion of the Hubbard model on the kagome lattice by the 
cluster dynamical mean-field theory, with the continuous- 
time auxiliary-field quantum Monte Carlo method as an 
impurity solver. We found that, in addition to charge 
and spin, the chirality degree of freedom becomes rele- 
vant in this system, and energy hierarchy appears among 
charge, spin, and chirality. Moreover, we showed that 
the quasiparticle state with strongly renormalized mass 
is stabilized at low temperatures, as a manifestation of 
the strong frustration. The heavy-quasiparticle state is 
formed by the entropy release associated with the emer- 
gent composite degree of freedom, the chirality. This 
chirality-driven mass enhancement provides a prototypi- 
cal example for the large mass behavior in frustrated sys- 
tems under strong correlation, and gives a clue to contro- 
versies on the heavy-fcrmion behavior of transition metal 
compounds such as LiV204. It is interesting to extend 
the study for more realistic models by considering the 
effect of carrier doping and orbital degrees of freedom. 
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